home *** CD-ROM | disk | FTP | other *** search
/ Magnum One / Magnum One (Mid-American Digital) (Disc Manufacturing).iso / d18 / nrpas13.arc / BESSK.PAS < prev    next >
Pascal/Delphi Source File  |  1991-05-01  |  407b  |  20 lines

  1. FUNCTION bessk(n: integer; x: real): real;
  2. VAR
  3.    tox,bkp,bkm,bk: real;
  4.    j: integer;
  5. BEGIN
  6.    IF (n < 2) THEN BEGIN
  7.       writeln('pause in routine BESSK');
  8.       writeln('index n less than 2'); readln
  9.    END;
  10.    tox := 2.0/x;
  11.    bkm := bessk0(x);
  12.    bk := bessk1(x);
  13.    FOR j := 1 TO n-1 DO BEGIN
  14.       bkp := bkm+j*tox*bk;
  15.       bkm := bk;
  16.       bk := bkp
  17.    END;
  18.    bessk := bk
  19. END;
  20.